DESY 04-158 



MS-TP-04-22 



^1" 

o 
o 

CO 



> 

o 
o 

On 
O 
^J- 
O 
-t— > 

Oh! 
<D 
43 



X 



Unquenched simulations with TV/ = 2 light quark flavours 
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The quark mass dependence of meson masses and decay constants in Nf = 2 QCD is studied at a fixed 
gauge coupling with different dynamical quark masses. Partially Quenched Wilson Chiral Perturbation Theory 
(PQWxPT) is applied to obtain the extrapolation to zero quark mass. It is shown, that in our analysis NNLO- 
terms play a more important role than C(a)-terms, which can be neglected. Also we compare our results with 
recent CP-PACS results. 



1. Introduction 

To have a reliable extrapolation from lattice- 
QCD results to physical parameter values, it is 
important to simulate (dynamical) quark masses, 
which are light enough in that sense that Chiral 
Perturbation Theory (%PT) P can be applied. 
If one enters this quark mass region, it is possi- 
ble to extract the Gasser-Leutwyler-constants [5] 
from lattice simulations, which then guide the ex- 
trapolation to physical quark masses. The chiral 
logarithms predicted by this low energy effective 
theory may serve as a sign that one reached suf- 
ficiently light quark masses. 

2. Results from PQWxPT 

Our simulations |3] were done at j3 = 5.1 with 
four different sea-quark (i.e. dynamical) masses. 
The lattice size is L 3 x T = i6 3 x 32. This work is 
a continuation of 0] were the lattice size was 16 4 
and only three different sea-quark masses were 
used. 

The lattice spacing was found to be a = 
0.195(4)fm, which gives a comfortably large vol- 
ume (L ~ 3fm). The dynamical quark masses 
(Wu,d) ar e in the range \m s < m u _d < \rn s , the 
pion masses are in 371McV < < 664MeV. 



*Talk presented by E. Scholz at Lattice 2004, Fermi Na- 
tional Accelerator Laboratory, June 21-26, 2004. 



The extraction of the Gasser-Leutwyler-constants 
is done by fitting ratios of masses and decay con- 
stants to the formulas predicted by PQWxPT pj]. 
We quote our results for the universal low energy 
scales, assuming fo = 93MeV: 

A 3 = f ■ 8.21(27) = 0.76(3)GeV, 
A 4 = f • 21.4(1.5) = 1.99(14)GeV. 

These fits also allow to extrapolate the pion decay 
constant to zero quark mass (for details cf. |3]): 

ZX l fo = 121(5)MeV. 

Comparing the new 16 3 x 32 data to the pre- 
vious 16 4 one, there are deviations up to 28% in 
the measured quantities like ro/a, pion and quark 
masses (which shows the largest deviations). For 
a more detailed discussion cf . [3] . It is remarkable 
that almost all these changes are substantially re- 
duced when comparing directly the ratios of me- 
son masses and decay constants as functions of 
ratios of the (PCAC-)quark masses, fig. n serves 
as an example. 

3. 0{o) vs. NNLO 

Since we are still working at a fixed gauge cou- 
pling /3, from our data we only have indirect in- 
formation about the discretization effects. That 
means we did not compare our results at differ- 
ent lattice spacings but fitted the 0(a)-effects in 
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Valence quark mass dependence of RRf, k=0.1770 
16 3 x32 ■ ■ 



Valence quark mass dependence of RRn+2RRf 



16 J x 32 lattice 

= 5.10 
k = 0.177 

Xs 3 0n2f= -0.0068(29) 
%=p/X s = 0.041(21) 




Figure 1. Simultaneous plot of 16 4 - and 16 3 x 32- Figure 2. Fit of NNLO- (solid) and (9(a)-formula 

data for RRf = /*f as a function of £ = (dashed) of RRn 4- 2RRf - 3. 
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the PQWxPT, which is possible due to the differ- 
ent functional dependence of these terms on the 
quark mass. Already in 0] we found the O(a)- 
contributions in the considered ratios and double 
ratios, parametrized by 

p 2W c sw 2B 

V = ~, P= 72 a, X s = ^ T rn q , 

AS Jo Jo 

to be small: rj < 0.1. When trying global fits 
(i.e. fitting the different ratios and double ratios 
of the valence mass dependence for all the four 
sea quark masses together) we still find the best 
fit at a small r\ « —0.01 (this would correspond 
to the effect of a quark mass of ks 0.5MeV), but 
a good fit with r\ = (no O(a)-terms) is also 
possible. To confirm this result, we looked at the 
combination of double ratios 

RRn + 2RRf - 3 = j . 

which vanishes in the continuum to NLO order 
and, therefore, it could be either described by 
NNLO- or C(a)-terms. The two fits are shown 
in fig. 0for our lightest sea quark mass. The x 2 
of the fit is 1.3 for NNLO and 7.2 for 0(a). This 
special example also suggests that the NNLO is 
more important in this fitting procedure. 



4. Comparison 

Other groups are also working (using different 
fcrmion actions, finer lattices etc.) on the same 
problem. For a review cf. [H], 0. The CP- 
PACS collaboration published results at almost 
the same lattice spacing (a = 0.2fm) with a set of 
light sea quark masses at 12 3 x 24 and 16 3 x 24 

First we compare the ratio ^"Jz^m VS- ( r o m Ps) 2 
(fig. |3J). Because we did not compute any renor- 
malization factor for the quark mass, we fitted 
our data to the CP-PACS data, i.e. assuming 
Z q = 1 for CP-PACS. The best fit gave for the 
qq+q data: Z q m 1.1142, which is shown in the 
plot. Of course the errors are large, but there is a 
reasonable agreement between the two data-sets. 

Next one may look at the sea-quark mass de- 
pendence of Rf = 4^-. By accident it hap- 

JR. Ft. 

pened that in the two data sets there is a point 
with almost exactly the same pion mass (romps): 
k = 0.14585 (CP-PACS) and k = 0.1765 (qq+q). 
This is used to set a common scale for compar- 
ing the sea-quark mass dependence of the pion 
decay constant, i.e. looking at the ratio Rf. In 
fig. 01 the data points are shown together with 
the (continuum formula) fit from The sec- 
ond (dashed) curve shows the fit, if one takes the 
value of A4 = 2.44(13)GeV from and uses the 
reference mass parameter \R = 35.8(3.3) p5J. The 
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CP-PACS, Z a =1 i Q ~ 
qq+q, £-1.1142 
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CP-PACS, 12,x24 ^ 
CP-PACS, 16>24 —A— 
qq+q, 16 3 x32 ' □ 
CP-PACS, continuum limit 
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Figure 3. Comparison of the ratio ^"Jzm vs ' 
(r mp S ) 2 between CP-PACS (0) and qq+q'(|3l)- 



qq+q fit gives A 4 = f .99(f4)GeV, which agrees 
within errors. The deviation of the CP-PACS 
data points to the continuum curves should be 
viewed as (large) C(a 2 )-effects. Furthermore the 
qq+q data-points show the right curvature due to 
the chiral logarithms, which is easily seen when 
connecting the two heaviest qq+q points by a 
straight line (dashed-pointed line in fig. 0J) . 

5. Discussion &: Outlook 

In the previous we showed the higher impor- 
tance of the NNLO-terms compared to 0(a)- 
corrections in our data and found a qualitative 
agreement with recent CP-PACS results. 

Recently the phase structure of lattice QCD 
has been examined by adding a twisted mass term 
([Hj and ref. therein, also [HJ], EU)- Besides the 
Aoki-phase scenario there exists a second scenario 
which results in a first order phase transition at 
K cHt . The possibility that our two lightest quark 
masses are already inside the metastable region 
can not be excluded, as also pointed out in [7j. 
This may be the reason for the observed strong 
non-linearity of the quark-mass dependence on 
(2k)~ 1 . We plan to investigate this question in 
the future. 

Apart from the study of the phase structure, 
the twisted mass formulation also facilitates nu- 
merical simulations at light quark masses. The 



Figure 4. Rf = £^ vs. a = ^ for CP-PACS 
and qq+q . Also shown are the (continuum 
limit) curves and a dashed-pointed straight line 
connecting the two heaviest qq+q points 



reason is the lower bound to the eigenvalue spec- 
turn provided by the additional mass term. This 
suggests to continue these types of simulations in 
the future with twisted quark masses. 
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